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The present contribution reviews recent advances made toward a micro- 
scopic understanding of superfluidity in nuclei using many-body methods 
based on the BCS ansatz and low-momentum inter-nucleon interactions, 
themselves based on chiral effective field theory and renormalization 
group techniques. 



1. Introduction 

1.1. Superfluidity in nuclei and BCS ansatz 

The structure of the nucleus and properties of extended nuclear systems 
strongly depend on their possible superfluid nature^ In nuclei, pairing 
impacts all low-energy properties of the system, e.g. masses, separation 
energies, radii as well as individual, rotational and vibrational excitation 
modes. The role of pairing correlations is emphasized close to the drip-lines 
due to the proximity of the Fermi surface to one- and two-particle emission 
thresholdsP In neutron stars, superfluidity also plays a key role, e.g. it 
impacts post-glitch timing observations^ or their cooling historyPl 

Despite the major role played by pairing in nuclei, its microscopic un- 
derstanding is rather poor. Given realistic inter-nucleon forces, questions of 
interest include (i) how much pairing is due to the direct attraction between 
nucleons on the one hand and how much it is due to the indirect coupling 
resulting from the exchange of medium fluctuations on the other, (ii) what 
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impact Coulomb and three-nucleon forces have. 

The key feature of any method based on the Bardeen-Cooper-Schriffcr 
(BCS) ansatz is to address explicitly the non-perturbative process associ- 
ated with the formation of quasi-bound (Cooper) pairs within the many- 
fermion system. This however done at the price of breaking the U(l) sym- 
metry associated with particle-number conservation, which eventually calls 
for its restoration. The latter is of particular relevance to mesoscopic sys- 
tems such as the atomic nucleus. In an ab-initio context, e.g., within the 
frame of many-body perturbation theory^ or self-consistent Green's func- 
tion theory^ the BCS rationale is incorporated by expanding the exact 
solution around an unperturbed state of the BCS or Bogoliubov type that 
already captures the key non-perturbative physics responsible for superflu- 
idity 



1.2. Low-momentum inter-nucleon interactions 

Establishing a realistic nuclear Hamiltonian, which is the basic precursor to 
any ab-initio many-body calculations, is a challenge for low-energy nuclear 
physics. The two-body sector has been intensively investigated and various 
interactions exist that reproduce nucleon-nucleon scattering phase shifts 
with xV-^dof ~ 1 in the elastic regime (up to about 300-350 MeV energy 
in the laboratory frame). The unsettled frontier is three- and higher-body 
forces.™ 

The development of chiral effective field theory (x-EFT) has made pos- 
sible^ to connect low-energy inter-nucleon forces to underlying Quantum 
Chromo Dynamics (QCD), whose relevant high-energy effects are renor- 
malized through fitted low-energy contact termsP^^ The main benefits of 
X-EFT are (i) to formulate the problem at hand in terms of relevant low- 
energy degrees of freedom (pions and nucleons) while retaining the (chiral) 
symmetry (breaking) of the underlying theory (QCD), (ii) to provide a 
consistent building of all relevant operators and (iii) to naturally explain 
the phenomenologically-observed hierarchy that makes two-nucleon (2N) 
interactions more important than three-nucleon (3N) interactions, which 
themselves dominate four-nucleon (4N) forces etc. Such a hierarchy relates 
to a power counting that organizes the infinite set of diagrams in the %-EFT 
LagrangiarfSl according to their scaling with (Q/A x ) u , where Q embodies 
typical low-momentum processes and degrees of freedom at play while A x , 
the so-called chiral-symmetry-breaking scale, denotes the hard scale char- 
acterizing omitted degrees of freedom and driving low-energy constants in 
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the Lagrangian. As such, %-EFT underlines that any nuclear Hamiltonian 
is effective and take the A x -dependent form 

H = T + V 2N + V 3N + ... (1) 

where V" 2N first contributes at leading order (LO) while V^ 3N only enters at 
next-to-next-to leading order (N 2 LO). 

Such nuclear Hamiltonians display several sources of non-perturbative 
behaviour that complicate many-body calculations. First are the strong 
virtual coupling between low- and high-momenturrE] modes driven by the 
non-observable short-range part of central and tensor forces. The second 
source of non-perturbative behaviour relates to the existence of weakly- 
and nearly-bound two-body states associated with 3 Si and x 5o partial- 
waves of the 2N interaction, respectively. Nearly-bound neutron-neutron 
and proton-proton states in the vacuum are the precursor of Cooper pairs 
that emerge at finite nuclear densities and that are central to the present 
discussion. 

Progress toward controlled calculations has long been hindered by the 
difficulty to solve the nuclear many-body problem in presence of the strong 
virtual coupling between low- and high-momentum modes. This has histor- 
ically been accepted as an unavoidable reality. Recently, EFT and renor- 
malization group (RG) methods' 14 * 15 ! have promoted a different view point 
based on the fact that the Hamiltonian (potential) is not an observable 
to be fixed from experiment. There rather exists an infinite number of 
Hamiltonians (potentials) capable of accurately describing the low-energy 
physicsPSI Starting from, e.g., a x-EFT Hamiltonian, one can take advan- 
tage of such a freedom to perform a (unitary) transformation U(A) with the 
aim of decoupling low-momentum modes from high-momentum ones. Doing 
so corresponds to keeping the physics (i.e. any true observable) invariant 
while transforming the Hamiltonian according to 

H(A) = U{A)HU f (A) =T + V W (A) + V m (A) + ... . (2) 

The RG transformation "lowers the resolution scale A" (< A x ) of the Hamil- 
tonian while preserving the original truncation error. Just as they depended 
on the original resolution scale A x , 2N, 3N,. . .AN interactions further de- 
pend on A while observables do not! 11 * 15 ^ The main benefit of the above 
transformation is that each component of H (A) typically becomes softer as 
A is lowered such that many-body calculations become more perturbative 
as far as the virtual coupling to high-momenta is concerned.^ One must 

a In the present context of %-EFT, "high- momenta" remain below A x . 
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however note that the source of non-perturbative physics associated with 
weakly- and nearly-bound two-body states, which is the focus of the present 
contribution, remains untouched by the lowering of A/ A. Even though 
Cooper pairs are tamed down as the density of the medium increases,^ 
they must be explicitly dealt with through non-perturbative techniques at 
and below nuclear saturation density. As discussed above, it is a key virtue 
of BCS-based methods to already account for the dominant effects of such 
a non-perturbative process through a simple zero-order ansatz. 

Other practical advantages of proceeding to a RG transformation over 
a reasonable range of A values relate to the fact that the unitarity of the 
transformation is not significantly compromised (and thus the physics not 
significantly altered) (i) by omitting induced many-body forces beyond a 
certain rank and (ii) by eventually truncating the size of the initial Hilbert 
space. Of course, the improved convergence pattern deriving from the lat- 
ter truncation and the preserved hierarchy of AN forces built into H must 
be thoroughly checked in actual many-body calculations! 20 " 23 ! As a mat- 
ter of fact, one can exploit the change of many-body observables with A 
to study the underlying physics scales and evaluate the incompleteness of 
approximate calculations or the impact of dropping many-body forces in 
the transformed Hamiltonian H(A). 

Two main classes of RG transformations are used to construct low- 
momentum interactions. For a schematic illustration of the evolution of 2N 
potentials' matrix elements in (relative) momentum space within each of 
this two methods, see Fig. 9 of RefP^ In the first approach, denoted as 
similarity renormalization group (SRG), the decoupling is achieved through 
a unitary transformation over the Hilbert (Fock) space defined in connec- 
tion with the original, e.g., %-EFT Hamiltonian. A momentum scale A 
that measures the extent of off-diagonal coupling is lowered through flow 
equations such that the potentials are driven toward a band diagonal form 
in momentum space. In the second approach, generically referred to as 
the View k approach, the decoupling is achieved by a transformation that is 
in fact not unitary over the original Hilbert (Fock) space. The renormal- 
ized potential is such that matrix elements beyond the lowered momentum 
cutoff A are set to zero. In practice, the low-energy part of the renormal- 
ized potentials are eventually very similar in both approaches such that the 
breaking of unitary in the V\ ow k approach has no influence on low-energy 
observables 
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1.3. Outline 

The present contribution is organized as follows. Section [5] discusses the ex- 
traction of the most basic observable related to pairing in nuclei. Section [3] 
reviews the recent use of realistic low- momentum interactions within the 
frame of the semi-empirical single-reference energy density functional (SR- 
EDF) method. Section|4]reports on the even more recent implementation of 
ab-initio self-consistent Gorkov Green's function calculations based on low- 
momentum interactions that constitutes a path towards a fully microscopic 
description of superfluidity in nuclei. 

2. Pairing information from the odd-even mass staggering 

The present contribution focuses on the most basic information related to 
pairing, i.e. the "pairing gap" extracted via three-point mass differences 

AW (AT) = L^i_ [i^+i _ 2E » + E"- 1 ] , (3) 

where N and Eq denote the number of nucleons and the ground-state 
binding energy, respectively. Extracting A( 3 )(7V) is motivated by the re- 
lation between the odd-even staggering of nuclear masses and the lack of 
binding of odd-even systems due to the presence of an unpaired nucleon.^ 
As such, A( 3 )(A^) is indeed dominated by the "pairing gap". Still, finite- 
difference formulas are contaminated by sub-leading contributions that are 
not related to the pairing gap 1 25 1 26 1 (see Eq. [13] below) . Consequently, and 
to avoid any mismatch when comparing theory to experiment, it is rec- 
ommended, whenever possible, to perform differences of actual theoretical 
masses rather than to use an approximate formula for the theoretical gap. 

3. Semi-empirical energy density functional calculations 
3.1. Elements of formalism 

The SR-EDF method implements the breaking of Z7 (1) symmetry and takes 
the form of an effective Hartree-Fock-Bogoliubov (HFB) formalism.^ The 
total energy is postulated under the form a functional Eq — £[p, k, k*} of 
the (symmetry breaking) one-body density p and pairing tensor n computed 
from an auxiliary product state |$o) of the Bogoliubov type, 

p ah = ($ |4c a |$o), Kab = (*o|Q,C a |$o) , (4) 
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where {c a } denotes an arbitrary basis of the one-body Hilbert space Hi. 
The minimization of £[p,K,n*] under the constrain that N — (&o\N\Q }, 
leads to solving effective Bogoliubov-De-Gennes equations^ 

fh-\ A 
^ -A* —h* + A 

where (U, V) M are the upper and lower components of Bogoliubov quasi- 
particle eigenstates whereas denotes the corresponding quasi-particle 
energies. The single-particle (h) and pairing (A) fields arc defined through 

tab + S^(cff) = t ab + ^bcl Pdc, (6a) 
cd 

Si 2 b (eff)^i^^ cd ^ , (6b) 

cd 

where t a b denotes the matrix elements of the kinetic energy operator. In 
Eq.[6j effective normal E n (eff) and anomalous E 12 (eff) self-energies, as well 
as effective particle-hole u p ' 1 and particle-particle v pp kernels, are introduced 
for interpretation purposes. The effective character of the EDF approach 
relates to the fact that £[p,n,K*] and E 9S (eff) are meant to re-sum cor- 
relations that go largely beyond Hartree, Fock and Bogoliubov diagrams 
calculated in terms of vacuum inter-nucleon interactions. 

Modern empirical parametrizations of existing, e.g. Skyrme or Gogny, 
EDFs provide a fair description of bulk and spectroscopic properties of 
known nucleiP^ However, they lack predictive power away from available 
data and a true spectroscopic quality, especially regarding the part that 
drives superfluidity. Consequently, efforts are currently made to empirically 
improve the analytical form and the fitting procedure of functionals, e.g. 
see RefsP2M22l f or recent attempts to pin down the isovector content of local 
pairing functionals. 

Along with improving the phenomenology at play, it is relevant to un- 
derstand the processes responsible for superfluidity in nuclei in microscopic 
terms, i.e. starting from vacuum inter-nucleon interactions. The lowest- 
order contribution to the two-particle irreducible pairing kernel v pp is pro- 
vided by vacuum interactions themselves, while higher-order terms include 
the induced interaction associated with the exchange of collective medium 
fluctuations between paired particles JHHH] A fundamental, yet unresolved, 
question relates to how much of the pairing gap in finite nuclei is accounted 
for at lowest order and how much is due to higher-order processes. 
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Awaiting for fully ab-initio calculations of mid-mass nuclei (see Sec.|4|), 
we have recently set up a semi-empirical SR-EDF scheme that combines 
empirical parametrizations of the single-particle field h with a pairing field 
A built at first order in 2]Np3S3 (including CoulomtP) and 3]NpE31 low- 
momentum interactions. The objective was to provide a partial answer to 
the fundamental question alluded to above. The next section summarizes 
the main results obtained following such an approach. 

3.2. Calculation scheme 

Our calculations start from the N 3 LO 2N potential (EM 500 MeV) of Ref.^B 
which is then RG-evolved using the Vj ow k procedure and a smooth n oxp = 4 
regulator^ down to A = 1.8 — 2.8 fm -1 . We complement rtP^ at each A 
by the leading N 2 LO 3N force, which is then averaged over the third nu- 
cleon following the procedure detailed in RefP^l Adding both contribu- 
tions, we obtain v pp at first order in vacuum interactions for various A 
values. To be consistent, the isoscalar and isovector effective masses of the 
Skyrme parametrization used to build v ph is constrained from Hartree-Fock 
calculations of neutron and symmetric nuclear matter based on the same 
low-momentum 2N and 3N interactions.^ Last but not least, odd-even nu- 
clei are computed through the self-consistent blocking procedure performed 
within the filling approximation^^!] such that A^(N) is indeed extracted 
through total mass differences. 

3.3. Results 

Figure Q] compares theoretical and experimental three-point mass differ- 
ences along several semi- magic isotopic/isotonic chains. Results obtained 
with and without 3N contributions to v pp are displayed. Using the 2N con- 
tribution only, neutron and proton pairing gaps are close to experimental 
ones for a large set of semi-magic spherical nuclei, although experiment is 
underestimated in the lightest systems. The addition of the first-order 3N 
contribution lowers pairing gaps systematically by about 30%. 

In Ref.p^l the A dependence of the pairing gaps was also studied as a 
way to estimate uncertainties due to short-range higher-order AN forces and 
to the incomplete many-body treatment. The analysis led to a theoretical 
uncertainty of order w 100 — 250 keV in all semi-magic chains investigated. 

All in all, the recent studies reported in RefsP^HlH have supported the 
following conclusions: (i) Coulomb repulsion is responsible for a signifi- 
cant proton anti-pairing effect and must be explicitly incorporated in the 
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Fig. 1. Theoretical and experimental A( 3 \N) along isotopic/isotonic chains based on 
low-momentum 2N and 3N interactions with A/A31MF = 1.8/2.0 fm" 1 . 



pairing kernel, (ii) it is essential to include the 3N contributions to the 
pairing kernel to reach a quantitative description of nuclear pairing gaps, 
(iii) first-order results leave about 30% room for contributions from higher 
orders, e.g. from the coupling of Cooper pairs to (collective) density, spin 
and isospin fluctuations, which is consistent with phenomenological cal- 
culations showing that induced interactions are overall attractive in nuclei 
[33j{36] ^j v ^ f u t ure investigations need to also compute the normal self-energy 
and higher-order contributions to the pairing kernel consistently from low- 
momentum 2N and 3N interactions. This is the objective of the next section 
to report on the first step in that direction. 



4. Ab-initio self-consistent Gorkov-Green's function theory 

The ultimate goal is to perform fully microscopic calculations of pairing 
properties in mid- and heavy-mass nuclei. One way to do so consists of per- 
forming ab-initio self-consistent Gorkov Green's function (SCGGF) calcu- 
lations^ based on low momentum interactions. While it is exact in the limit 
where self-energy diagrams are summed to all orders, SCGGF reduces to the 
(ab-initio) HFB approximation at lowest order. In practice, calculations are 
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the result of a compromise, i.e. one must implement a tractable truncation 
scheme that approximates the exact solution well enough. In that respect, 
switching from conventional hard-core potentials to low-momentum inter- 
actions is instrumental as it makes finite-order schemes qualitatively and 
quantitatively viableP^ 

Recently, V. Soma, C. Barbieri and myself have implemented SCGGF 
theory at second order on the basis of low-momentum 2N interactions ) 48 l 49 l 
These are the first ever ab-initio calculations of their kind in finite nuclei. 
They constitute a first step towards a quantitative investigation of induced 
interaction effects from a purely ab-initio perspective. To be of quantitative 
interest, such calculations need to be extended to 3N interactions and to 
more advanced truncation schemes allowing for the coupling of the Cooper 
pair to collective fluctuations. These constitute our mid-term objectives. 

4.1. Elements of formalism 

Let basis {c\ } splits into two blocks that can be mapped onto each other 
through time-reversal. Let us then introduce a partner basis {c],} through 

<t(t)=vA(t), Ca(t) = VaCa(t) , (7) 

which corresponds to exchanging the state a by its time-reversal partner a 
up to the phase rj a . By convention a = a and r\ a r/a = — 1 . 

In Gorkov formalism, one targets the ground state |\l/o) of the grand- 
canonical-like potential fl = H — [iN, where fi is the chemical potential and 
N the particle-number operator, having the number N = (\E , o|-^|^'o) of 
particles on average. In order to access the complete one-body information 
contained in |^o)i one must introduce a set of four Green's functions, known 
as Gorkov propagators.^ Defining an "annihilation" column vector and a 
"creation" row vector through 

C a (f) = , Ct(t) = ( c t( t ) c a (t)) , (8) 

one can write the four propagators in matrix representatiorPHl through 

iG ab {t,t') = {^ \T\c a (t)Cl(t')\\^ } = i (9) 

V GH(t,t>) Gf b (t,t>) J 

Self-consistent, i.e. dressed, propagators are solution of Gorkov's equation 

/ T + £ n (u;)-A E 12 (w) 
^ S 21 (w) -T + £ 22 H + A 
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whose outputs are Gorkov's amplitudes (U, V)^ and corresponding pole en- 
ergies uifj,, in terms of which the four propagators can be expressed.^ Equa- 
tion [TU] generalizes Eq. [5] in the sense that normal and £ 22 (w)) and 
anomalous (£ 12 (w) and £ 21 (w)) irreducible self-energies act here as energy- 
dependent potentials. At first order in vacuum interactions, Eg. 1101 reduces 
to an ab-initio^! HFB equation with normal and anomalous self-energies 
accounting for Hartree-Fock and Bogoliubov diagrams, respectively. Pro- 
ceeding to an actual calculation relates to truncating the diagrammatic 
expansion of the self-energies T," (w). As opposed to perturbation the- 
ory, the expansion involves skeleton diagrams expressed in terms of dressed 
propagators solution of Eq. (fTU|) . This key feature of self-consistent Green's 
function methods allows the re-summation of self-energy insertions to all 
orders and makes the method intrinsically non-perturbative and iterative. 
Eventually, the total energy is computed via the Koltun-Galitskii sum rule^ 

E ° = E f / duj [tab + 00 Sab] ' 

where the integration is performed over a closed contour in the upper half 
of the complex plane. Extensive details regarding both the formalism and 
the computational scheme can be found in RefsPSHSl 

Computing odd-even nuclei requires in the present context to perform 
Gorkov calculations for a state |^o) having an odd number-parity quan- 
tum number. This is however beyond the scope of the present work. The 
next best approximation consists in keeping an even number-parity state 
while accounting for the self-consistent blocking of a quasi-particle within 
the filling approximation. 4 ^ Such an approximation remains however to be 
formulated within the frame of Gorkov-Green's function formalism. Conse- 
quently, the next best estimate of the ground-state energy of an odd-even 
system is obtained througrP^SI 

E»mE»*+ U », (12) 

where Eq * is the energy of the odd nucleus computed as it were an even 
one, i.e. as a fully paired even number-parity state having an odd number 
of particles on average, while uip denotes the lowest pole energy extracted 
from Eq.Hnifor that calculation. For even N one simply has Eq = Eq *, i.e. 
Eq * provides the energy curve on which both even and odd nuclei would 

b This is at variance with the effective character of Eq. \E\ in which energy-independent 
fields h and A are meant to effectively account for correlations that go beyond strict 
Hartree, Fock and Bogoliubov diagrams. 
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lie in the absence of odd-even mass staggering. With such a decomposition 
of the energy, Eq. |3] becomes 

A(3) (7V) «L2L^_ + AF(iV) . (13 ) 

The second derivative of E^ 1 * is smooth with N but provides A( 3 \N) 
with a rapidly oscillating contribution^ because of the factor (— 1) N that 
comes with it in Eq. (|13|) . Clearly, such a contribution is not related to 
the pairing gap. The second contribution to A^ 3 \N) relates specifically to 
the unpaired character of the odd nucleon and does extract, in open-shell 
nuclei, the pairing gap at the Fermi energy^^ 

a f^s \ for N odd ,„ 

«- 1 + V 1 )/2 for N even ' < 14 > 

4.2. Calculation scheme 

We performed self-consistent second-order calculations i.e. first- and 
second-order diagrams are included in the computation of the self-energies 
(see Figs. 1-3 of Ref. 48 ), of Ca isotopes. Such a truncation scheme con- 
stitutes a Kadanoff-Baym ^-derivable approximation, which automatically 
ensures the exact fulfilment of conservation lawsP^ The x _ EFT 2N poten- 
tial is RG-evolved down to A = 2.1 fm _1 using the Viowfe procedure and 
a sharp regulator. In addition to the direct (i.e. first order) contribution 
to the pairing kernel incorporated (along with the corresponding diagram 
from the 3N interaction) in the semi-empirical EDF calculations reported 
on in Sec. [3l the present calculation includes the coupling of the Cooper 
pair to non- collective density, spin and isospin fluctuations. Such a test case 
calculation is performed on a single-processor and uses a restricted model 
space of 7 major harmonic oscillator shells. As such, results presented here 
must only be taken as indicative. Converged calculations require a multi- 
processors architecture and will be reported on in Ref: 49 

4.3. Results 

Figure [2] compares theoretical and experimental three-point mass differ- 
ences from 37 Ca to 51 Ca. Results obtained at first and second-order are 
displayed. We observe three main features (1) the second-order contribu- 
tion to normal and anomalous self-energies generates a slight decrease of 

c See Figs.Q] 
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A( 3 )(7V), (2) gaps at first and second order account for about half of ex- 
perimental values, (3) the staggering of A^(A^) is inverted compared to 
experiment. All such features will be discussed in details in RefP^ on the 
basis of better converged calculations performed using larger single-particle 
model spaces. Let us however make a few tentative comments regarding 
each of these three points. 

(1) The coupling to non-collective fluctuations seems to have little influence 
on A( 3 '(vV) and to have the tendency to slightly suppress it. Although 
one should wait for converged calculations and for a thorough analysis 
of the result before drawing any conclusion, the fact that the collective 
character of the fluctuations the Cooper pair couples plays a key role 
is expectedP^HiH Treating collective fluctuations within the present 
ab-initio setting is a challenge that will be addressed in the mid-term 
future. 

(2) The omission of the 3N force in the computation of the normal self- 
energy results in a significantly too small effective mass. This prevents 
the pairing kernel from fully expressing its strength and eventually 
generates too low gaps. The second-order contribution to the normal 
self-energy from the 2N interaction does not compensate (enough) for 
this defect such that incorporating the 3N force contribution is manda- 
tory to obtain a realistic effective mass and meaningful predictions of 
the pairing gapfl It particular, having a realistic effective mass will 
feedback on the second-order contribution to the anomalous self-energy 
discussed in point (1) in a way that remains to be seen. Note that the 
indirect enhancement of the gaps associated with the 3N force contribu- 
tion to the normal self-energy will counteract its repulsive contribution 
to the pairing kernel discussed in Sec. 13.31 

(3) The inverted staggering of A^ 3 ) (N) emphasizes (see first term of Eq.[T3j) 
the wrong curvature of Eq * obtained in the present calculation as one 
adds neutrons.^ Such a feature relates to the lack of saturation of nu- 
clear matter, and correspondingly to a wrong asymmetry energy, gener- 
ated by low-scale 2N interactions when omitting the corresponding 3N 
interaction.^ Within each shell, the energy per particle increases as one 
adds neutrons whereas it should decrease. Increasing the scale A of the 
2N interaction partly corrects for such a wrong pattern. Still, resolving 
this issue satisfactorily necessitates, independently of the value of A, 



d In the context of EDF calculations discussed in Sec. [3] the realistic effective mass is 
obtained empirically by fitting appropriate free parameters of the Skyrme EDF. 
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N 

Fig. 2. Neutron three-point mass difference in Ca isotopes. Theoretical values are from 
first and second-order SCGGF calculations based on a low-momentum 2N interaction 
with A = 2.1 fm — . See text for details. 

the explicit treatment of 3N interactional within the present ab-initio 
setting. This again constitutes one of our two mid-term objectives. 



5. Conclusions 

The present contribution has reviewed recent advances made toward a mi- 
croscopic understanding of superfluidity in nuclei using many-body methods 
based on the BCS ansatz and low-momentum inter-nucleon interactions. 
Steps towards a truly ab-initio treatment of pairing in medium-mass nuclei 
have been briefly highlighted. 
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